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(Anomalous) X-ray Pulsars 

V. M. Kaspi^*, F. P. GavriiP 

'^Physics Department, McGill University, Rutherford Physics Building, 3600 University St. 
Montreal, QC, Canada H3A 2T8 

We review the observational properties of the class of young neutron stars known as "anomalous X-ray pulsars," 
emphasizing the tremendous progress that has been made in recent years, and explain why these objects, like 
the "soft gamma repeaters," are today thought to be young, isolated, ultrahigh magnetic field neutron stars, or 
"magnetars." 



1. INTRODUCTION 

Prior to the commissioning of BeppoSAX and 
the Rossi X-ray Timing Observatory in 1996, 
the so-called "Anomalous" X-ray Pulsars (AXPs) 
were considered very mysterious sources, because 
the energy source for their bright X-ray emission 
was unknown. At the time, there were only 3 
known members of this class. They were distin- 
guished by having periods in the narrow range 
6-9 s, showing approximately steady spin-down, 
and having softer spectra in general that those 
seen in accreting X-ray pulsars. All were known 
to lie within 1° of the Galactic Plane, and inter- 
estingly, one source, IE 2259-1-586, was known 
to reside in the supernova remnant CTB 109. 
AXPs as a class were identified as having mod- 
est X-ray luminosities, in the range ^ 10'^^ — 



The leading model to explain 



the AXPs was that they were accreting neu- 
tron stars, though with properties very differ- 
ent from the bulk of established accreting X- 
ray pulsars, including the absence of any evi- 
dence of a companion I* va n Paradiis et all Il994 
[Mereghc tti & Stella, 1995). 

The situation post-BeppoSAX and especially in 
the latter years of RXTE is very different and 
much clearer. The basic phenomenology of the 
sources is now well mapped out. Here, we sys- 
tematically review the most important proper- 
ties of this class of objects, which now includes 
5 and possibly 8 sources (see Tables [T] and EJ, 



*Steacie Fellow, Canada Research Chair 



and summarize why today, accretion models are 
strongly disfavored. Rather, the magnetar model, 
in which AXPs are isolated young neutron stars 
powered by a decaying ultrahigh magnetic field, 
provides the most compelling explanation for the 
unusual AXP source properties, as it does for 
an equally as exotic class, the soft gamma re- 
peaters (S GRs). AXPs have also been reviewed 
recently bv lMereghetti etal 1 ll2002l) . and magne- 
tars in general have been reviewed recently by 
.KaspL (,2003a-bl . 

2. TIMING PROPERTIES OF AXPS 

Since their discovery, AXPs have been known 
to be spinning down. Unlike most known ac- 
creting X-ray pulsars, no evidence was seen, in 
nearly two decades of timing, for any extended 
spin-up. However, some deviations from sim- 
ple spin-down were observed. AXP 1E2259-I-586 
showed a handful o f possible very short li ved 
spin-up events (e.g. iBaykal & Swank', '19961) as 
did IE 1048-5937 (e.g. Oosterbrock ct al., 1991). 
These were noted by various authors and were 
suggested to be due to accretion torque vari- 
ations (e.g. fea vkal fc Swank, 1996), ghtches 
llHevl & Hernau stL Il999() and magnetar radia- 
tive precession ijMelatosl fl997|) . However, with 
sparse observations consisting of a frequency mea- 
surement every few years and rarely more often, 
determining the origin of the apparent deviations 
from simple spin-down could not be done. 

In order to address this problem, a program of 



1 



Table 1: Spin parameters for AXPs. 
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IE 1048-5937 


> 2.7 




6.45 


~ 3.81 


~ 5.0 


~ 55 


~ 2.7 


2 


RXS 1708-4009 


~ 8 




11.00 


1.86 


4.6 


5.4 


9.4 


3 


IE 1841 045 


5.7-8.5 


Kes 73 


11.77 


4.16 


7.1 


9.9 


4.5 


4 


IE 2259+586 


3 


CTB 109 


6.98 


0.0483 


0.59 


0.55 


230 


5 


AX J1845.0-0258* 


~ 8 


Kes 75 


6.97 










6 


CXOU J0110043.1-721134* 


57 




8.02 










7 


XTE J1810-197* 


~ 10 




5.54 


1.15 


2.6 


26 


7.6 


8 



(*) not confirmed; (f ) sec Ozcl, Psaltis & Kaspi 2001 for a discussion on distance estimates for the confirmed AXPs; References: (1) 
Gavriil & Kaspi 2002; (2) Kaspi et al. 2001; (3) Kaspi & Gavriil 2003; (4) Gotthelf et al. 2002; (5) Woods et al. 2003; (6) Torii et 
al. 1998; (7) Lamb et al. 2003; (8) Ibrahim et al. 2003. 



Table 2: Spectral parameters for AXPs. 



Source 


riH 


r 


kT 


Lx 


/pi {%V 


Ref. 




(xl022 cm-2) 




(keV) 


(erg s"^) 






4U 0142+61 


0.88 


3.3 


0.42 


3.3 X 10^4 


~ 88 


1 


IE 1048 5937 


1.0 


2.9 


0.63 


3.4 X 1034 


~ 80 


2 


RXS 1708-4009 


1.49 


3.11 


0.45 


6.8 X 10^5 


~ 73 


3 


IE 1841-045 


2.0 


2.26 




2.3 X 10^^ 


100 


3 


IE 2259+586 


0.93 


3.6 


0.41 


1 X 10^^ 


~ 50 


4 


AX J1845.0 0258* 


9.0 


4.6 




7.4 X 1034 


100 


5 


CXOU J0110043.1-721134* 


0.14 




0.41 


1.5 X 10^5 





6 


XTE J1810-197* 


1.05 


3.75 


0.668 


1.6 X 10^6 


~ 70 
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(*) not confirmed; (j) contribution of the power-law component to the total flux, see Perna et al. 2001 for further discussion; 
References: (1) Juett ct al. 2002; (2) Tiengo et al. 2002; (3) Mereghetti et al. 2002; (4) Patel et al 2001; (5) Torii et al. 1998; (6) 
Lamb et al. 2003; (7) Ibrahim et al. 2003. 
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regular monitoring of the 5 confirmed AXPs by 
RXTE was initiated in 1998. The goal was to 
accomplish phase-coherent timing, in which ev- 
ery rotation of the neutron star is counted on 
time scales of months to years. Phase-coherent 
timing is done regularly for radio pulsars, and is 
effective with any periodic source in which the 
periodicity is very stable, or at least changes rel- 
atively slowly. This turned out t o apply nicely 
to the AXPs teasp i et all Il999l) . For exam- 
ple, the RMS phase residual for IE 2259-f586 
in ^5 yr of timing (pre- June 2002) is under 
2% of the pulse period, following the removal 
of a model having on ly three free parameters 
i Gavriil fc Kas^ 120021 hereafter GK02). Phase- 
coherent timing on long time scales has now 
been accomp l ished for AXPs RXS J1708-4009 
i Kasr^i et a,l.l . Il999t). 4U 0142-^61 (GK02) and 



IE 1841-045 l|Gotthelf et al.Ll2002^ and indicates 
these sources are capable of great rotational sta- 
bility. This stability argues against an accretion 
origin of the X-rays, since most accreting sources 
show much li i gher l evels of torque noise (but see 
iBavkal et all l200ll) . The stability is compara- 
ble in some cases (particularly 4U 0142-1-61 and 
IE 2259-f 586) to that seen in young radio pulsars. 
Together with the much noi s ier tim ing properties 
of SGRs fe.g. IWoods et al.Lll999l) . this provides 
support for a continuum of timing noise proper- 
ties in the radio pulsar, AXP and SGR popula- 
tions, in line with the magnetar model. 

However, one AXP, IE 1048-5937, is a much 
noisier rotator than the others, so much so 
that phase-coherent timing cannot be accom- 

Jilishcd over more than a few months at a time 
Kaspi et al.L l200l(l . More detailed observations 
of the source reveal that its spin-down rate can 
change on time scales of weeks, and by large fac- 
tors (see Fig. Gavriil & Kaspi, in preparation) . 
This behavior is reminiscent of that se en in SGRs 
1806-20 and 1900-^14 ( Woods et al.l l2Q02). 

Thus, deviations from simple spin-down in 
AXPs appears to come in three flavours: (i) 
glitches and subsequent recovery; (ii) low-level 
stochastic variations having a "red" spectrum, 
similar to the "timing noise" seen in radio pul- 
sars; and (iii) large, short-time-scale variations 
which preclude phase connection. The origin of 



the latter two in particular is unknown. The 
low-level variations in radio pulsars may be re- 
lated to crustal superfluid effects such as "mini- 
glitches" , or may, in some cases, result from long- 
term recoveries from glitches that p receded the 
commencement of the observations. lArras et al.l 
l|2003l) have recently suggested that the larger- 
scale torque variations arise from angular momen- 
tum transfer from a superfluid core. Such a core, 
they argue, also results in a reduction in the inte- 
rior temperature that could make the crust more 
brittle, hence result in greater burst activity as 
seen in the SGRs (and possibly IE 1048-5937; 
see 
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Figure 1. Top: Long-term frequency history of 
IE 1048-5937 (after Kaspi et al. 2001). The 
heavy lines represent intervals over which phase- 
coherent timing was possible. Bottom: Same 
data but with the linear trend removed. 



2.1. Glitches 

Because phase-coherent timing counts every ro- 
tation, it determines spin parameters with high 
precision. This permits sensitivity to glitches 
having fractional amplitudes as low as ^ 10~^. 
The first AXP glitch was detected in RXS 
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J1708-4009 te aspj et all |200(^), and had frac- 
tional amplitude 6 x 10~^, and an increase in 
the magnitude of the spin-down rate of ^1%. 
These glitch properties are similar to those seen in 
Vela-like radio pulsars. Interest ingly, this source 
glitched again ^1.5 yr lat er ijKaspi fc Gavriil 
1200.4 iDairOsso et all Eool . However, the sec- 
ond glitch was much larger, with fractional fre- 
quency change 4 x 10^®, and a significant post- 
glitch recovery in which nearly all of the glitch 
relaxed on a time scale of '^50 days. The fre- 
quency history of this source is shown in Figure|21 
Neither glitch was accompanied by any obvious 
radiative changes. 




Modified Julian Day 

Figure 2. Frequency history of RXS J1708-40G9 
showing the two very different glitches. The 
top plot shows the overall frequency evolution, 
while the bottom plot shows the same but with 
the long-term spin-down trend removed, as well 
as measured frequencies. The best-fit model is 
based on a phase-coherent analysis (after Kaspi 
& Gavriil 2003). 



The second discovered AXP glitch was in 
IE 2 259-H586 l|KasT)i et all l200,4 IWoods et all 
l200.1l . Unlike the RXS J1708-4009 glitches, this 
one occurred simultaneously with (or possibly a 
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Figure 3. Frequency history of IE 2259-f586 
around the time of its 2002 outburst based on a 
phase-coherent analysis. The top plot shows the 
frequency evolution around the glitch, along with 
measured frequences. The middle panel shows 
frequency residuals, while the bottom shows 
phase residuals (after Woods et al. 2003). 



few hours before - see Woods et al. 2003) a major 
outburst in which over 80 X-ray bursts were de- 
tected in just a few hours, in addition to sudden 
order-of-magnitude increases in the pulsed and 
unpulsed flux, significant pulse profile changes, 
and an infrared enhancement (all discussed be- 
low) . This represents the first neutron-star glitch 
ever observed to be accompanied by significant 
radiative changes, and clearly indicates a major 
event that simultaneously affected both the inter- 
nal and external structure of the star. Roughly 
20% of the glitch recovered on a time scale of 
weeks, and in doing so resulted in the stellar spin- 
down being a factor of > 2 greater than its pre- 
outburst value (Fig. O. This is unprecedented 
in radio pulsars, and suggests that just following 
the glitch, the neutron star superfluid was actu- 
ally spinning slower than the crust, with the ob- 
served subsequent enhanced spin-down a result of 
angular moment transfer from the crust back to 
the superfluid (Woods et al. 2003). Additionally, 
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there was evidence that the ghtch may have been 
resolved in time, on a time scale of ~2 weeks. 

Glitches are definitely expected in t he ma gne- 
tar model fe.g.lThompson fc DuncanL Il996|) . As 
pointed out by iKaspi e t al. (20001), at least in 
principle, an accreting source can undergo a spin- 
up glitch since the latter results from an inter- 
nal angular momentum transfer from superfluid 
to crust regardless of the nature of the external 
spin-down torque. However, one would not ex- 
pect simultaneous bursts in an accretion scenario, 
as one might in the magnetar model. 

The observed AXP glitches provide a very plau- 
sible explanation for the historically observed 
spin-down deviations at least in IE 2259-1-586 
(Baykal & Swank 1996). The picture is still not 
clear for IE 1048-5937, however. 

3. AXP X-RAY PULSE PROFILES AND 
PULSED FRACTIONS 

AXP pulse profiles are, like those of the SGRs, 
broad, with large (^80%) duty cycles, and gen- 
erally significant harmonic content (e.g. GK02). 
The profiles show energy dependences that vary 
from source to source. A possible trend of greater 
energy dependence for profiles with higher har- 
monic content was identified by GK02, who also 
showed that in general, AXP pulse profiles are 
very stable. 

However, in 2002 June, simultaneously with the 
detection of the glitch and X-ray bursts, the pulse 
profile of IE 2259-1-586 underwent significant 
changes, on time scales from hours to days (Kaspi 
et al. 2003; Woods et al. 2003). The profile had 
relaxed back to its pre-outburst morphology by 
^2 we eks following the outburst. Ilwasawa et alJ 
ill992l) observed an apparent change in the X-ray 
pulse profile of IE 2259-1-586, in which the rel- 
ative amplitude of the two peaks in the profile 
changed between observations made in 1989 and 
1990. This can be explained as being due to an 
outburst having occurred just before the 1990 ob- 
servation (see [jSJ. 

AXP pulsed fractions vary from source 
to source, with the highest being ^0.8 for 
IE 1048-5937 and the lowest being -0.1 for 
4U 0142-1-61. Some, but not all, are energy depen- 
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Figure 4. Pulse profile changes in IE 2259+586 
seen around the time of the 2002 June outburst 
(after Woods et al. 2003). 



dent, and those which are vary differently with 
energy. For a summary of AXP pulse d fractions 
and t heir energy dependences, see lOzel et al.l 
ll200lh . It is not clear whether the pulsed frac- 
tions are time variable in general. However, the 
pulsed fraction of IE 2259+586 clearly changed at 
the time of its 2002 outburst: immediately post- 
outburst, the pulsed fraction decreased from its 
quiescent level of ^0.23 to ^0.15, however it re- 
covered fully after 3 days (Woods et al. 2003). 

4. X-RAY SPECTRA 

X-ray spectra of AXPs generally require two 
components. These are usually taken to be a 
thermal blackbody component with a power-law 
tail. The measured spectral parameters of the 
known AXPs are given in Table |21 The spectra 
as a class are softer than those of the SGRs in 
quiescence. The softest source in that class is 
SGR 0525—66; its spectral parameters are actu- 
ally softer than those of IE 104 8-5937. which 
among other things, prompted *Kulka rni et al.l 
(2003) and Kaspi et al. (2001) to suggest these 
sources may be transition objects between the two 
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classes. 

In the context of the magnetar model, the spec- 
tra can be understood as follows. The ther- 
mal component is emerging from the stellar sur- 
face, a result of heating of the interior by ac- 
tive magnetic field decay dThompson fc DuncaiA 
119951: iThompson fc DuncanL Il99fil) . The ther- 
mal spectrum is thought to deviate significantly 
from that of a blackbody, because of the effects 
of the stellar atmosphere, as well as the large 
magnetic field, which results in different opaci- 
ties for different photon polariz ations, as well as 
on QED vacuum polarization ('Ho & Lai', '2004 
2001; Zanc ct al., 2001; Ho & Lai, 200^ 



Ozel l2003|) . The thermal spectrum is hardened 



relative to a blackbody of the same temperature 
due to the non-grey atmosphere, although vac- 
uum polarization counteracts this slightly. As ob- 
servers fit the thermal component with a black- 
body, some portion of the non-thermal compo- 
nent may result from the atmospheric distor- 
tion. However, this portion is probably small. A 
more promising origin of the non-thermal emis- 
sion is external resonant Compton scattering of 
thermal seed photons by magnetospheric currents 
(jrhomDson et al., 2002). 

The X-ray spectra were, pie- Chandra and 
XMM-Newton, hoped to hold direct evidence 
for the high magnetic field v ia features such 
as electron cyclotron lines (e.g. IHo fc Lail I2OOI1 
IZane et all 1200 it lOzel I2003HHo fc Lai l2003l) 
Of course an electron cyclotron line in a ~ 10^^ G 
field might look similar to a proton cyclotron 
line in a ~ 10^^ G field. In any case, no such 
lines have been seen in spite of some high spec- 
tral, and i n some cases , temporal reso lution ob- 
servations llPatel et a l.'. '2001*; 'juett et 
iTiengQ et al.ll2002HWoods et al...200: 



;t all 



5. AXP X-RAY FLUX VARIABILITY 

Strong flux variability pie-BeppoSax and 
RXTE was reported for IE 2259+586 and 
IE 1048-5937 (Baykal & Swank 1996; Ooster- 
broek et al. 1998). Flux variations of a fac- 
tor of 5-10 were reported, albeit from different 
instruments, having different spectral responses, 
with some imaging and some not. However Iwa- 



sawa et al. (1992) reported a brightening of a 
factor of r^2 in a 1990 GINGA observation of 
IE 2259-1-586 compared with an observation in 
1989. They noted that the 1990 pulse profile was 
also significantly different than that observed pre- 
viously, with different relative peak amplitudes, 
and different peak shapes. Furthermore, the mea- 
sured 1990 spin period was fractionally shorter 
by ^ 3 X 10^^ compared with what the previous 
spin-down rate would have predicted. 

In ~5 yr of monitoring using the PGA on 
RXTE, GK02 found no evidence for such flux 
variations in any AXP. This was consistent with 
what was found by Tiengo et al. (2002) in a 
comparison of past observations of IE 1048—5937 
with recent XMM-Newton data. The overall re- 
cent lack of variability in AXPs thus appeared 
discrepant with the historical record. 

The 2002 June outburst of IE 2259-1-586 ap- 
pears to have solved this conundrum, at least 
for this source. Simultaneous with the bursting 
were increases in the and persistent fluxes by a 
factor of > 10 (Kaspi et al. 2003; Woods et 
al. 2003), which mostly decayed on a time scale 
of days, but which has left an X-ray afterglow 
in which the pulsed flux is still a factor of ~2 
greater than the pre-outburst value a year since 
the outburst (Fig. . The total energy in excess 
pulsed and persistent emission during the short- 
decay-time-scale enhancement was 3 x 10^^ erg 
(2-10 keV), while that in the extended afterglow 
is much more, 2 x 10''^ erg (Woods et al. 2003). 

The rapidly decaying flux enhancement seen in 
IE 2259-1-586 could be due to a transient sur- 
face hot spot. During the rapid initial flux de- 
cay, the blackbody radius was smaller than at 
all other times, the temperature was higher, and 
the pulse profile was clearly different, support- 
ing this picture. Alternatively, it could have 
been magnetospheric, as a large current density 
will be excited in the magnetosphere above re- 
gions of strong crustal shear. The short-lived af- 
terglows detected after intermediate SGR bursts 
have a simple explanation as the cooling of a 
pair-r ich surface lay e r hea ted by a high-energy 
flare l|Ibrahim et However, no such 

flare was seen for IE 2259-1-586. This is prob- 
lematic also for explaining the long-time-scale af- 
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Figure 5. The pulsed flux history of IE 2259+586 
(2—10 keV) around the time of the outburst, in- 
dicated by the sharp spike (after Woods et al. 
2003). 



terglow. In SGRs, bulk heating of the crust can 
power an excess heat flux from its surface for 
a year or more, and has been proposed as the 
explanation for the quasi-power-law flux decay 
seen in SGR 1900 +14 (.Lvubarsky et al. 2002 ') 
and SGR 1627-41 llKouvehotou et all l2003i) "ln 
each case, an initial deposition of 10'*'' ergs was 
assumed, consistent with the detection of an ini- 
tial giant soft gamma-ray flare; this was unseen 
for IE 2259+586. For a more detailed discussion 
of the possible origins of the enhanced emission, 
see Woods et al. (2003). 

We note that the combination of the observed 
flux enhancement, glitch and pulse profile change 
in IE 2259+586 observed by Iwasawa et al. 
(1992) using GINGA are consistent with an out- 
burst similar to that observed in 2002 June having 
occurred days/ weeks prior to their 1990 observa- 
tion (see Sjnil- This offers an estimate of a crude 
burst rate of two every '^20 yr. 

This observed variability associated with out- 
bursts also makes the two transient AXP 
candidates (see Table AX J1845-0258 

i Vasisht GottheHlh 9971 ). and XTE J1810-197 
Ibrahim et al.[ 120031) easier to understand. 
These two objects have both shown factor of > 10 
increases in their fluxes. For XTE J1810-197, the 



flux decreased slowly after its appearance, in con- 
cert with its spin-down rate, not unlike the behav- 
ior seen in IE 2259+586 post-outburst (Ibrahim 
et al. 2003; Woods et al. 2003). Such tran- 
sient AXPs suggest a large population of quies- 
cent AXPs exists in the Galaxy. 

6. X-RAY BURSTS 

The first discovery of bursts from AXPs came 
from the RXTE/PCA monitoring observations of 
IE 1048-5937. Two faint bursts, separated by 
^2 weeks, w ere detected in ^42 5 ks of exposure 
over ~5 yr 1 Gavriil et al.L 120021) . These bursts 
very much resemble SGR bursts. Specifically, 
their fast rise times, short durations, hard spec- 
tra relative to the quiescent emission, fluence and 
probably clustering, are all SGR burst hallmarks. 
The origin of the bursts could not unambiguously 
be proven to be the AXP, given the large PCA 
field-of-view, and the absence of any other radia- 
tive or spin change in the source. Intriguingly, 
the first burst's spectrum was not well fit by a 
continuum model, showing evidence for a strong 
emission line at ~14 keV. 

Not long after the reporting of the above 
two bursts, a major outburst consisting of 
over 80 bursts was detected from the direction 
of IE 2259+586 fortuitously during a regular 
RXTE/PCA monitoring observation in 2002 June 
(Kaspi et al. 2003). These burst s were very simi- 
lar to those of SGRs ijGavriil et al .. 2003). Specif- 
ically, like the SGRs, the AXP burst durations 
follow a log-normal distribution which peaks at 
99 ms, the differential burst fluence distribution 
is well described by a power law of index —1.7, 
the burst fluences are positively correlated with 
the burst durations, the distribution of waiting 
times is well described by a log-normal distribu- 
tion of mean 47 s, and the bursts are generally 
asymmetric with faster rise than fall times. 

However, there were some notable differences 
between the AXP and SGR bursts that may be 
clues to the physical differences between the two 
source classes (Gavriil et al. 2003). Specifi- 
cally, the AXP bursts exhibit a wider range of 
durations and, unlike SGR bursts, occur prefer- 
entially near pulse maxima; the correlation be- 
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tween burst fluence and duration seen for SGRs 
is flatter than for SGRs; the AXP bursts are on 
average less energetic than are SGR bursts; and 
the more energetic AXP bursts have the hard- 
est spectra - the opposite of what is seen for 
SGRs (Gavriil et al. 2003). Furthermore, in stark 
contrast to SGRs, the energy detected in bursts 
(6 X 10^^ erg, 2-60 keV) was much smaher than 
that in the post-outburst persistent flux enhance- 
ment (2 X 10^-^ erg, 2-10 keV). This could indicate 
bursting activity that was missed by our observa- 
tions and the gamma-ray monitors, although the 
latter would have easily detected SGR-like bursts 
having the missing energy (Woods et al. 2003). 
This "quiet" outburst strongly suggests there are 
many more such objects in the Galaxy than was 
previously thought, as is also indicated by the 
transient AXP candidates (see 

Overall, the properties of the outburst in 
IE 2259-1-586 argue that the star suffered a ma- 
jor event that was extended in time and had two 
components, one tightly localized on the surface 
of the star (i.e. a fracture or a series of fractures) 
and the second more broadly distributed (possi- 
bly involving a smoother plastic change). The 
glitch points toward a disturbance within the su- 
perfluid interior while the extended flux enhance- 
ment and pulse proflle change suggest an excita- 
tion of magnetospheric currents and crustal heat- 
ing. The very rich data set provided by this out- 
burst should be very useful in constraining phys- 
ical properties of the affected neutron star. 

7. OPTIC AL/IR OBSERVATIONS AND 
VARIABILITY 

Of the five confirmed AXPs, four now have 
secure or possible optical/IR counterparts. The 
first optical/IR detection of an AXP w as made of 
4U 0142-K61, bv lHnlleman et a,lJ tOOCt) . They ar- 
gued that the source, which had i? ~ 25 mag, was 
too dim to b e from an accretion di sk. This was 
confirmed by iKern fc Martini l)2002() who showed 
that this source is pulsing with the X-ray pe- 
riod and a 27% pulsed f raction, much tog high 
to be reprocessed light. iH^iUeman et al.l lj200ll) 
discovered a possible near-IR {Kg = 21.7 mag) 
counterpart to IE 2259+586. This was con- 



firmed by Kaspi et al. (2003) who found the 
source to have brighted by a factor of ~3 three 
days after its 2002 outburst, but faded by a fac- 
tor of ~2 1 week later. This source appears 
to have stayed brighter than pre-ou t burst, how- 
ever, t w o months after the outburst lllsrael et al.l 
l2003b>) . IWang & Chakrabartvl l| 2002l) reporteTa 
likely near-IR counterpart to IE 1048—5937 hav- 
ing K = 19.4 mag, and dete cted in multiple wave- 
bands. 'Isra el et al.l l)2002|) reported significant 
IR variability in the source on a time scale of 
^50 days. They speculated the variability might 
be related to bursting activity from this source 
(Gavriil et al. 2002), in an alogy with the vari- 
ability seen in IE 2259-1-586. Ilsrael et all l)2003a|) 
reported the detection of a possible IR counter- 
part to RXS J1708-4009. 

These detections are interesting for two rea- 
sons. First, as noted above, in two cases vari- 
ability is observed; for IE 22594-586, it is likely 
that it is associated with its 2002 June outburst. 
Such variability may therefore prove to be an im- 
portant observational model constraint, although 
currently its origin in the magnetar model is un- 
clear. Most likely the optical/IR emission is a 
product (or, given the large X-ray to optical/IR 
luminosity ratios, L^/Li^ > 500, a byprod- 
uct) of radiation processes in the outer magneto- 
sphere, and therefore is sensitive to the changes 
in the current structure induced by magnetic re- 
configurations. Continued monitoring for corre- 
lated optical/IR and torque variations seems war- 
ranted, especially as a way of testing the pro- 
posed "twisted magnetosphere" model of mag- 
netars (Thompson, Lyutikov & Kulkarni 2002). 
Second, as pointed out bv Ilsrael et alJ |2003^, 
the spectral energy distributions of the optical/IR 
and X-ray emission show that the former is much 
too bright to be the simple extrapolation of the 
blackbody component of the X-ray spectrum. 
However, it is fainter than the extrapolation of 
the X-ray power-law spectral component, so call- 
ing it an "excess" may be premature. 
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8. CONCLUSIONS AND OPEN QUES- 
TIONS 

Since the 1996 commissioning of BeppoSAX 
and RXTE, our overall picture of AXPs has 
changed dramatically. The number of likely 
AXPs has nearly tripled, and our understanding 
of these unusual sources' properties has improved 
tremendously. Perhaps the single most impor- 
tant discovery is that the apparent resemblance of 
AXPs with SGRs noted by Thompson & Duncan 
in 1995 is more than skin deep: with the discovery 
of bursts from AXPs, the two source classes are 
now united unambiguously. Our next challenge is 
to learn how to extract physically interesting in- 
formation from AXP and SGR observations. In 
this sense, their study is still in its infancy, with 
observations ahead of theory. 

This said, there are obvious possibilities for 
fruitful observational investigation of AXPs. 
First, it would be nice to have a direct mea- 
surement of the inferred high magnetic field. As 
no X-ray spectral features have been forthcom- 
ing, another avenue is needed. X-ray polarization 
observations are an excellent possibility, partic- 
ularly for the brightest AXPs. The detection of 
such polarization, in addition to confirming the 
high magnetic field, would be the first demon- 
stration of the birefringence of the vacuum, as 
predicted by QED. In the shorter term, glitches in 
AXPs may offer a practical method of constrain- 
ing the structure and physics of these objects. 
The simulteneity of the IE 2259-f 586 glitch with 
its outburst and associated radiative changes, 
we suspect, is telling us a lot about the stellar 
structure. Continued patient timing of these ob- 
jects has the potential to reveal correlations be- 
tween glitch properties like amplitude and relax- 
ation time scales with radiative properties, which 
will help us understand properties of the highly 
magnetized crust and superfluid interior. Opti- 
cal/IR observations also offer hope of constrain- 
ing magnetar outer magnetosphere processes, al- 
though its origin is not yet clear. Finally, there 
is the open question of the radio pulsar/ AXP 
connection. Recently, several radio pulsars hav- 
ing inferred magnetic fields higher than that of 
IE 2259-1-586 have been discovered, yet with no 



evidence for any AXP-like X-ray em ission (e.g. 
^ Pivovaroff et al., McLaiighlin et al., .20m . 

This is puzzling. It may simply reflect the fact 
that the magnetic field measured by P and P is 
approximate only, in which case the discovery of 
more such radio pulsars should eventually result 
in the identification of the "missing link." 

Funding for this work comes from NSERC 
(via a Discovery Grant and Steacie Supplement), 
the Canada Research Chair Program, NATEQ 
(via Team and Observatoire de Mont Megantique 
Grants), CIAR (via a Fellowship), and NASA 
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